Magneto-optical properties of Europium hexaboride 
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Full potential spin-polarized self-consistent electronic structure calculations have been performed 
on LaBe, EuBe and YbBg using Linear MufBn-Tin Orbital Method. The electronic structure of 
EuBe shows a small overlap near the Fermi level. This is in disagreement with the recent ARPES 
study. The reflectivity spectrum shows a plasma resonance at about 0.3 eV in EuBe. Absorption 
starts at about 1.5 eV and it is due to 4/ to 5d transition. Both these features are reflected in 
the magneto-optical Kerr spectrum. At 0.3 eV, there is a rotation of +8 to -10 and at 1.5 eV, a 
rotation of about 0.5° occurs. The reasonable agreement with MOKE experiment implies that the 
inclusion of Drude correction is absolutely essential for the treatment of EuBg. 

PACS numbers: 



I. INTRODUCTION 

Divalent hexaborides have been at the focus of atten- 
tion for the last couple of years because of their anoma- 
lous electronic and magnetic properties. Divalent rare 
earth hexaborides, EuBg and YbBg behave quite differ- 
ently due to the half and full filled / orbital of Eu and 
Yb respectively. YbBg is nonmagnetic whereas EuBg or- 
ders fcrromagnctically ' at 15.1 K accompanied by a large 
negative magneto resistance and a significant blue shift" 
of the refiectivity plasma edge.' A second phase transi- 
tion occurs at 12.7 K with a broad peak in the specific 
heat or an anomaly in the resistivity' Neutron diffrac- 
tion experiment indicates that magnetization increases 
very sluggishly with decreasing temperature in the ab- 
sence of magnetic field. ^ The spectacular discovery of 
ferromagnetism in alkaline earth hexaborides. La-doped 
CaBg, CaBg, SrBg and Ba-doped CaBg have aroused 
deep interest in these compounds.''''''' Curie temperature 
of alkaline hexaborides is above 600 K, much higher than 
EuBg. The introduction of chemical pressure by substi- 
tuting Ca for Eu in EuBg lowers ferromagnetic transition 
temperature to 5.5 K. 

The intriguing magnetic and transport properties of 
hexaborides arise from the complex electronic proper- 
ties. Spectroscopic measurements indicate that EuBg 
may be a narrow gap semiconductor, in unanimity with 
early" band structure calculation. However the resis- 
tivity measurements revealed metallic behavior below 
300 K.[l()] Shubnikov-de Haas, de Haas- van Alphen^^ 
and refiectivity^" measurements suggest that EuBg is a 
semimetal with a very low carrier concentration whereas 
LaBg is a semimetal with a higher carrier concentra- 
tion. Model Hamiltonian calculations' ' claim that elec- 
tronic properties of EuBg may be explained by consid- 
ering the fact that it is a semimetal. The anomaly in 
electrical resistivity suggests that it is close to metal- 
insulator transition. However the recent optical experi- 
ments observe that there is a sizable increase in charge- 
carrier concentration and/or a reduction of the effec- 
tive mass of the itinerant charge carriers as it is cooled 
through the ferromagnetic transition temperature. Hall 



effect measurements indicate that with the increase in 
magnetic order, the effective carrier concentration in- 
creases significantly. ' ' Angle-resolved photoemission and 
resonant inelastic X-ray scattering experiments reveal 
that EuBg is a semiconductor with a gap of 0.8 eV.[15] 
Magneto-optical Kerr measurements suggest a Lorentz- 
Drude model in order to account for a large resonance in 
its Kerr rotation spectrum. '^^' There is stress upon the im- 
portance of the interplay between the charge carriers and 
the localized-electron interband transitions. All these ex- 
periments imply that EuBg has a metallic or a semimctal- 
lic ground state. Hence the ground state of EuBg is not 
only controversial but an open problem. 

First principle band structure calculations based on 
local density approximation (LDA) yield semimetal with 
small overlap between the valence and conduction band 
at the X point of Brillouin zone.^''^'"' Again these re- 
sults are within the limitations of density functional the- 
ory (DFT). The failure of LDA in describing electronic 
structure of semiconductors is well known. LDA is also 
proving to be inadequate for some s and p systems and 
grossly wrong for the strongly correlated systems like 
the transition metal and rare earth (RE) compounds. 
Thus corrections to LDA such as Self-interaction^'* and 
inclusion of on-site Coulomb interaction ( orbital depen- 
dent potential ) are required to obtain the exact elec- 
tronic structure of correlated systems. The band struc- 
ture calculation performed within weighted density and 
GW approximations-" has successfully produced a gap in 
CaBg. Strong on-site Coulomb interactions among the 
4/ electrons in rare earth systems describe their optical 
properties satisfactorily.^^'^^'^"^'^^ Very recently LDA -I- U 
calculations- ' have been performed on EuBg which con- 
cludes that its ground state is half metallic semimetal. 
Kunes et. al. discussed the origin of magnetism in terms 
of exchange interaction between conduction electrons and 
local 4/ magnetic moments. There have been as many 
advocates of the semiconducting nature of Eu+^Bg as 
of its semimetallic ground state. FLAPW calculations 
have been reported for EuBg without considering the 
strong correlations and optical properties of EuBg are 
not calculated.^' Magneto-optical spectra provide use- 



ful information about interplay between eleetronic strue- 
ture and the magnetization of systems. The optieal and 
magneto-optical properties of hexaborides have not yet 
been studied from the band structure calculation. We 
have investigated systematically optical and magneto- 
optical properties of ReBg (Re = La, Eu and Yb). The 
main purpose of our present work is to observe the effects 
of /" occupation of Re, where n=0, 7, 14, on optical and 
magneto-optical spectra and to compare with the exper- 
imental results. '^'''"'''^' 



II. CALCULATIONAL DETAILS 

Linearized band structure calculations are fairly ac- 
curate within the limitations of various approximations 
used in bypassing the many-body problem. Full poten- 
tial (FP) self-consistent spin-polarized band structure 
calculations have been performed on rare earth hexa- 
borides ReBg by the Linear MufSn-Tin Orbital Method 
(LMTO)-'' using the FP-LMTO code developed by S. Y. 
Savrasov et. al.."" Since La has empty 4/ shell, the cal- 
culations are performed within simple LSDA (local spin 
density approximation). The DFT is modified in order 
to include the strong correlations among the / electrons 
in EuBg and YbBg. In the LDA + U method '", the LDA 
energy functional is modified by removing the LDA / — / 
interactions and adding the strong on-site Coulomb inter- 
actions among the / electrons. The main idea of the LDA 
-f U method used in this calculation is that the LDA gives 
a good approximation for the average Coulomb energy of 
/ — / interactions E^u as a function of the total num- 
ber of / electrons. Subtracting this average energy from 
the LDA total-energy functional and adding orbital- and 
spin-dependent contributions yield the exact total energy 
within the mean field approximation. The detailed de- 
scription of the method is given in Ref . [24] . 



TABLE L Table contains a, the lattice parameter, muffin- 
tin radius of Re ion, Smt, structural parameter x, Coulomb 



parameter U and 


Density of States at Fermi level, N{Ef) 


a 
(a.u.) 


Smt 
(a.u.) 


X 


U 
eV 


NiEp) 
(St./eV) 


LaBe 7.920 
EuBe 7.855 
YbBe 7.832 


3.96 
3.93 
3.92 


0.202 
0.202 
0.202 


0.0 
9.0 
9.0 


0.86 
0.36 
0.0 



The value of the parameter U = (i^°) used in the 
present calculation is 9 eV for both the divalent hexa- 
borides in order to place the occupied / states at about 
2 eV below Fermi level, i?F-[-^)l] The other Slater integrals 
F^, F^ and F^ for / electrons have been used from the ex- 
perimental paper by Thole et. al. '- to calculate the ma- 
trices Umm' and Jmm' ■ The skeletal structure of these cu- 
bic hexaborides may be viewed as a CsCl type lattice with 
RE ion occupying the Cs sites and CI replaced by Bq oc- 
tahedral cages. Octahedral cages of Boron ions arc linked 



together along the three Cartesian directions. The struc- 
ture is formally known as CaBg type crystal structure and 
has a space group Pm'im{Oh). An internal parameter x 
determines the ratio between inter and intra octahedron 
B-B distances. Structural parameters x and lattice pa- 
rameter a are given in Table I. Table I also provides the 
parameter U and density of states at Ep. In SrBg, a lin- 
ear relationship between x and energy gap is found. For 
X > 0.206, band gap assumes a positive value. ^' . How- 
ever we use X = 0.202, as reported^' for EuBg. The 
exchange-correlation potential in the LDA was calcu- 
lated with Vosko-Wilk-Nussair parametrization. The ba- 
sis consists of 6s,5p, Sd, 4/ states of cation and Boron: 
{n)s and {n)p, where n refers to the principal quantum 
number. Charge density, density of states and the mo- 
mentum matrix elements were calculated on a grid of 242 
k points in irreducible Brillouin zone and the k space in- 
tegration was performed using the Tetrahedron method. 
We have treated orbital momentum quantum number, 
mi as a parameter. Eu ion is considered to be in its diva- 
lent state and calculations were performed with the 4/^ 
configuration. 



III. RESULTS AND DISCUSSION 

The self-consistent full potential spin-polarized band 
structure of EuBg is shown in fig. 1. 




FIG. 1: Spin-polarized band Structure of EuBs. 

Typical small overlap between the 5d conduction band 
and the Boron p valence bands is seen at the X(IOO)- 
point. The lower Hubbard 4/ band is situated in be- 
tween 2 and 3 eV below the Fermi level, Ep- The upper 
hubbard band lies about 5.5 eV above Ep. The bonding 
2p orbitals of Boron mainly constitute the valence band 
near Ep. The bottom of the conduction band also has 
some Boron 2p character. Hence the bottom of the con- 
duction band and the top of the valence band is formed 
from the hybridized Eu 5d and Boron 2p states. This pic- 
ture is fairly consistent with the XPS and BIS results ''' 
and previous FLAPW calculation- ' . It was suggested by 



Kiincs ct. al. that a ferromagnetic GW calculation may 
result in a gap opening up in spin-up channel only. How- 
ever this picture is contrary to most experimental results. 
On comparison with band structure of LaBg, as given in 
fig. 2, one finds that the Boron 2p bands lie much below 
Ef- The Fermi level has Eu f — d hybridized bands lo- 
cated on it. Apparently, this proves that LaBg has much 
more free carriers than EuBg. The unoccupied 5d band is 
at about 4 eV above Ep, which is about 1 eV lower than 
that in EuBg. The 4/ band lies 2 eV above Ep in LaBg. 
The band structure of YbBg reveals a zero gap, with zero 
overlap at X-point. The 5d states move further up in en- 
ergy in the conduction band. All the 14 electrons in the 
4/ state are occupied and the occupied 4/ band has a 
width of 2 eV. The position of the occupied 4/ band is 1 
eV deeper than that obtained from experiments. " The 
blue shift' ' of the 5d band is evident from the movement 
of the 5d conduction band towards higher energy from 
LaBg to YbBg. 
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FIG. 3: Band Structure of YbBe. 




FIG. 2: Band Structure of LaB 
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The basic difference between the divalent hexaborides 
and LaBg is that the overlap at the X-point is larger in 
LaBg and the overlap is between La d — f hybridized 
and Boron 2p states. Here it is to be kept in mind that 
the calculation for LaBg is performed within LSDA. The 
LSDA band structure is in agreement with previous'*'' 
band structure calculation. In the divalent hexaborides, 
there is small or no overlap between the Eu (or Yb) 
5d and Boron 2p states. Band structure calculation in 
EuBg [ I 7] reveal that the electronic structure is sensitive 
to the dimensions of the Boron octahedra. Choosing a 
larger [x > 0.206) value of the internal parameter x opens 
up a small band gap. It is apparent that even this large 
value is insufhcient in producing the gap of about 0.8 
eV.[L^] There is no available band structure calculation 
on YbBg as far as we know. 

Figure 4 depicts the angular-momentum decomposed 
density of states (DOS). The region near Ep comprises 



of Boron p and Eu d states in EuBg . The fundamental 
electronic structure of YbBg is very similar to that of 
EuBg since both the cations are divalent. The occupied 
/ state split into the 4/^/^ and 4/^/^ states in YbBg. 
The former lies at about 3.5 eV and the latter at about 
2 eV below Ep . The Boron p states lie deeper in LaBg 
than in the divalent hexaborides. 
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4: /-projected Density of States of ReBg (Re = La, Eu, 



The Fermi level is predominantly occupied by the La 
5d and 4/ electrons. It has no partially filled / orbital. 
However 5d states are partially occupied, as is evident 
from fig. 4. That the carrier concentration is very small 



in the divalent hexaborides is also evident from Table I. 
The decreasing value of DOS at Ep from LaBg to YbBg 
bears testimony to it. However one cannot compare the 
two divalent hexaborides directly from this data. The 
position of the 4/ band is fairly consistent with the XPS 
and BIS results in the divalent hexaborides. 

As mentioned earlier, LaBg and YbBg are nonmag- 
netic. The spin-polarized calculations yield zero total 
moment for both of them. EuBg has a total magnetic 
moment of 6.98 ^b- Orbital polarization is negligible for 
both the ions. Major contribution to the total moment 
is from the spin polarization of the 4/ electrons. Eu+^ 
has f"^ configuration. It has completely filled 4/ spin-up 
states and completely empty spin-down states. Hence 
the spin moment should approach 7 hb- Orbital mo- 
ment should be quenched for an / configuration. Our 
calculational results are consistent with the atomic ana- 
log for Eu. The value of /_* from neutron diffraction 
experiments '" is 7.3 ± 0.5 fiB- The large error bar comes 
mainly from uncertainty in the absorption. But the total 
moment from the present calculation is within that error 
bar. 

Optical conductivity tensor has been calculated from 
the standard relations given by Reim et. al.. ''^ For accu- 
racy calculations have been performed upto 30 eV. The 
imaginary component of diagonal optical conductivity 
i'^xx) ''^^'^ the real component of off-diagonal conductiv- 
ity {crl ) have been calculated using the Kramcr's-Kronig 
transformation. A description of the formalism is given 
in Ref. no. [24]. Optical and magneto-optical spectra 
have been calculated using the LMTO band structure 
and further, a Lorentzian broadening of about 0.0068 eV 
has been used in order to obtain the best agreement with 
the experimental results in the divalent Re hexaborides. 

Reflectivity spectra for these compounds are shown in 
figure 5. The inclusion of Drude correction reproduces 
the experimental reflectivity spectrum" of LaBg. An 
absorption at low energy occurs due to Eu 5d to Eu 4/ 
transition. Due to this transition, a small disagreement 
occurs at low energy. This can be removed if the 4/ states 
are pushed up by 1 eV. Addition of a constant potential 
to the Hamiltonian''" which acts only on the 4/ states 
of La has been successful in reproducing the reflectivity 
spectrum in La monochalcogenides. A large dip in re- 
flectivity or the plasma edge is successfully reproduced 
near 2 eV. However for the divalent hexaborides, a very 
small value of plasma frequency in the Drude correction 
is required to reach an agreement. The small Drude cor- 
rection required; howsoever small it may be, strengthens 
the claim that the divalent Re hexaborides are semimet- 
als with low carrier concentration. Due to smaller num- 
ber of free carriers, plasma edge effect occurs around 0.3 
eV and 0.2 eV in Eu and Yb hexaborides respectively. 

The calculated and experimental optical conductiv- 
ity spectra are depicted in figure 6. Drude correction 
brings in greater agreement with experiment as far as 
the plasma minimum is concerned. It occurs at about 2 
eV in LaBg. But the minima is hardly distinguishable 
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FIG. 5: Calculated reflectivity spectra of ReBe and their cor- 
responding experimentaI(Ref. 1,8, ■5.")) reflectivity spectra. 



for divalent hexaborides. The broadness of the minimum 
emphasizes their low carrier density. The first structure 
in EuBg at 2 eV is due to 4/ to 5d transition. This small 
absorption is missing in the spectrum for LaBg and hence 
the claim in the last line is justified. This structure shifts 
very slightly to higher energy in YbBg since the unoccu- 
pied 5d states shift to higher energy in the conduction 
band. The numerous structures till 6 eV are due to tran- 
sition from Boron 2p state to Eu 5d state. From fig. 4, 
it is obvious that the largest peak in EuBg will occur at 
7 eV and it will be due to transition from the Boron 2p 
main peak in valence band to the Eu bd main peak. This 
peak shifts to lower energy in LaBg and to higher energy 
in YbBg due to similar blue shift of the 5d state in the 
conduction band. Around 8 eV, 5d to 4/ transition and 
2p to 5d transition structures overlap. Beyond 9 eV, the 
structures arise solely due to transitions from the Boron 
2p states to the 5d peaks of Eu. 

The electron energy loss spectrum, as given in figure 7, 
supports the fact given above. The first peak at 0.2 eV 
for YbBg and that at 0.3 eV for EuBg, are surely due to 
the plasmon of the conduction electrons. The first peak 
is found at 1.8 eV in LaBg. This reduction in the energy 
position is due to the fall in the carrier concentration in 
divalent hexaborides. 

Magneto-optical Kerr effect is a very complex phe- 
nomenon. Here a very subtle interplay between spin- 
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FIG. 6: Calculated real component of diagonal optical con- 
ductivity ReBe and their corresponding experimental (Ref. 
35) counterpart. 



FIG. 7: Calculated energy loss spectra and the corresponding 
experimental(Ref. 3")) spectra for the ReBe. 



orbit coupling and spin polarization leads to MO signals. 
However, in most cases, as in semimetallic magnetic rare- 
earth compounds, giant signals are obtained due to the 
plasmon resonance. The most important prescription for 
obtaining large MO response is a large plasmon resonance 
followed by a strong interband transition. Inspite of large 
SO interaction and a fair plasma resonance, we fail to see 
any kind of MO response in YbBg shown in the left panel 
of figure 8, as expected for non-magnetic materials. This 
is also in agreement with experiment."'' The same hap- 
pens for LaBg. This proves that magnetism plays a very 
crucial role in MOKE. EuBg, with a total moment of 7 
/is, is totally spin polarized. The calculated spectrum, 
as given in figure 8, shows a jump from -1-8° to —10*^ 
at 0.3 eV and followed by a small signal of about 0.5° 
at 1.5 eV. The former is due to large plasma resonance 
and the latter is due to 4/ to SO split hd band tran- 
sition. On comparison with the experimental spectrum, 
the first signal finds agreement. The second one is slightly 
suppressed and occurs at 0.5 eV higher in energy. This 
degree of agreement implies that the Drude correction 
plays a very important role. Obviously this strengthens 
the fact that EuBg is semimetallic. Moreover, magnetism 
plays a crucial role along with SO coupling. 

Magneto-optical^'' response in undoped Eu hexaboride 
has been found to be giant at high fields and low tem- 
perature. It has been pointed out by Broderick et. al. 
that Drude component is essential in reproducing the 
overall features of the experimental Kerr rotation spec- 



trum. The resonance at 0.3 eV is caused by the response 
of the itinerant charge carriers to magnetism while the 
Kerr response at 1 eV is associated with the magneto- 
optical response of the 4/-5(i interband transitions. On 
doping EuBg with isoelectronic Ca, the sharp onset of 
plasma edge is broadened and the related Kerr resonance 
is wiped out. This can be due to the large damping of the 
itinerant charge carriers, which is in turn is due to disor- 
der in the cation sublattice.^'' The intensity of the struc- 
ture at 1 eV also diminishes with Ca doping since the 
contribution of the 4/ electrons of Eu decreases. These 
experimental observations are, hence, consistent with the 
electronic structure calculation of EuBg, presented here. 



IV. CONCLUSION 

Self-consistent spin-polarized band structure calcula- 
tions are presented for LaBg, EuBg and YbBg. They 
reveal the metallic nature of the first and the low car- 
rier density nature of the divalent ones. Small overlap 
at the X-point implies semimetallic nature of EuBg and 
YbBg. This is in disagreement with the recent ARPES 
measurements in EuBg. Present calculations lead to rea- 
sonable agreement with optical conductivity, reflectivity 
and magneto-optical Kerr spectrum. The successful re- 
production of the MOKE spectrum in EuBg implies that 
magnetism, spin-orbit coupling and large plasma reso- 
nance are responsible for the large magneto-optical re- 
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FIG. 8: Calculated Kerr spectra of LaBe, EuBe and YbBe and the experimental Kerr rotation of EuBs at 1.6 K and at 1 T 
field. (Ref. 27) Calculated Kerr spectra for LaBe and YbBe are magnified by a factor of 2 and 4 respectively. 



sponsc in this system. The MO response at about 1.5 cV 
is due to 4/ to 5d transition. The inclusion of strong on- 
site Coulomb repulsion is absolutely essential for its re- 
production. Most importantly, Drude-like nature at low 
energy is a very strong evidence for the semi-metallic na- 
ture of EuBg and this is strongly supported by the present 
calculations. 
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